Nonlinear Evolution of q = 1 Triple Tearing Modes in a Tokamak Plasma 
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In magnetic configurations with two or three q = 1 (with q being the safety factor) resonant 
surfaces in a tokamak plasma, resistive magnetohydrodynamic modes with poloidal mode numbers 
m much larger than 1 are found to be linearly unstable. It is found that these high-m double or 
triple tearing modes significantly enhance through nonlinear interactions, the growth of the m = 1 
mode. This may account for the sudden onset of the internal resistive kink, i.e., the fast sawtooth 
trigger. Based on the subsequent reconnection dynamics that can proceed without formation of 
the m — 1 islands, it is proposed that high-m triple tearing modes are a possible mechanism for 
precursor-free partial collapses during sawtooth oscillations. 
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PACS numbers: 52.55.Fa, 52.35.-g, 52.55.Tn 

The profile of the safety factor q(r) (which measures 
the magnetic field line pitch) contains information about 
the instability characteristics of magnetically confined 
plasmas in toroidal or helical systems with respect to 
current-driven magnetohydrodynamic (MHD) instabili- 
ties. In particular, in the presence of magnetic surfaces 
where q = 1, a mixing of the plasma inside these surfaces 
was observed. This instability is thought to be closely 
related to internal disruptions, generally known as saw- 
tooth oscillations, which strongly affect the quality of en- 
ergy and particle confinement. In view of the desired ap- 
plication to thermonuclear fusion reactors such as ITER, 
a detailed understanding of these internal large-scale in- 
stabilities is necessary. 

A heuristic model proposed by Kadomtsev |l| success- 
fully explains overall phenomena associated with a full 
sawtooth crash. In this model, a perturbation with helic- 
ity h = m/n = 1 (m being the poloidal and n the toroidal 
Fourier mode number), which is in resonance with the 
closed field lines on the q = 1 surface, quenches the hot 
core region inside the q = 1 surface through magnetic 
reconnection. This kind of relaxation, generally known 
as the m — I internal resistive kink instability (in short, 
the m = 1 mode) has been observed to exhibit an abrupt 
onset, which is called a "fast trigger." However, a satis- 
factory explanation of this fast trigger is not yet known. 
Another unresolved problem is the possibility of a par- 
tial sawtooth collapse, where the m = 1 mode saturates 
before the reconnection of the core is completed, so only 
an annular (off-axis) region undergoes mixing |2(. 

During the evolution of a tokamak plasma subject to 
sawtooth relaxation oscillations, multiple q — 1 resonant 



surfaces may arise temporarily jjj]. When this occurs in 
configurations with a hollow current profile (qo > 1 with 
go being the q at the magnetic axis) o = 1 double tearing 
modes (DTMs) can become unstable 0,|a,|fj whereas, for 
a centrally peaked current profile (go < 1), q = 1 triple 
tearing modes (TTMs) may arise |f|. 

In this letter, it is demonstrated that some of the 
phenomena associated with the sawtooth oscillations in 
tokamaks can be explained by the nonlinear evolution of 
DTMs or TTMs. Using nonlinear numerical simulations, 
it is shown that, in the presence of multiple q — 1 res- 
onant surfaces, rapidly growing high-m DTMs or TTMs 
can enhance the growth of the m = 1 mode and later 
generate electromagnetic turbulence in the annular re- 
gion surrounded by the q = 1 resonant surfaces. Based 
on these observations, it is shown that the fast trigger of a 
sawtooth crash as well as precursor-free partial collapses 
during sawtooth relaxation oscillations can be accounted 
for by the nonlinear evolution of TTMs. 

The set of equations we use is the reduced magnetohy- 
drodynamic (RMHD) equation in the zero-beta limit in 
a cylindrical geometry. In normalized form the RMHD 
model can be written as 
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where V is the magnetic flux function, <f> the stream 
function (electrostatic potential), j = — V^f/ 1 the ax- 
ial current density and u — Vj_0 the vorticity, essen- 
tially following the standard notation (cf., e.g., Ref. 0). 
The time t is normalized by the poloidal Alfven time 
t Hp (time scale for dynamics in an ideal magnetized 
plasma) and the radial coordinate by the minor radius a 
of the plasma column. The resistivity profile is given by 
fj(r) = j(r — 0,t — 0)/j(r,t — 0). As to the Lundquist 
number Sn p = tr/th p [where tr = /io/(?7oa 2 ) is the re- 
sistive time scale and rjo the resistivity at r = 0] and 
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FIG. 1: A safety factor profile q{r) unstable to TTMs. Verti- 
cal lines (dashed) indicate the locations of the q = 1 resonant 
surfaces, r s j < r S 2 < r B s, and q — 1 is indicated by a horizon- 
tal (dotted) line. 
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FIG. 2: Linear growth rate spectra of TTMs for the q profile 
given in Fig.^for the Lundquist numbers Sh p = 10 6 , 10 7 and 
10 9 and constant Prandtl number Pr = Sh p v = 1. 

viscosity v (normalized by o 2 /th p ), ^Hp = 1/1 = 10 is 
used, unless stated otherwise. The constant source term 
Eq compensates the resistive dissipation of the equilib- 
rium current. Using a quasi-spectral code with a finite- 
difference radial mesh, single helicity (h = m/n = 1, 
< to < 127) nonlinear simulations of Eq. |Q and @ 
were performed with an initial condition corresponding 
to a flow-less equilibrium with perturbed flux function 
(if>m>o(t = 0) = 10 , m — n). An ideally conduct- 
ing wall was taken as the boundary condition. Although 
the RMHD model encompasses only part of the phys- 
ical effects involved in sawtooth oscillations (cf., e.g., 
Refs. |2>l3)IS|)i it suffices for the present purpose of fol- 
lowing the fundamental dynamics of DTMs and TTMs. 

The important new feature of linear instability that 
is addressed here is the fact that configurations with 
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FIG. 3: Linear growth rate spectra of DTMs for q profiles 
with two q = 1 resonant surfaces located a distance D12 apart, 
obtained with Sn P = v~ x = 10 6 . In the three cases shown, 
D12 = 0.05, 0.1 and 0.2, respectively. 
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FIG. 4: Radial structure of linearly unstable TTM eigen- 
modes for (a) the m = 1 modes of ip, (b) the m = 13 modes of 
ip, (c) the m = 1 modes of cfi, and (d) the m = 13 modes of <f> 
obtained for Sh p = v~ X = 10 6 . For a given m, the eigenmode 
extending to the resonant radius r s , is denoted by Mi (to). 
Vertical lines (dotted) indicate resonant radii. 



multiple q — 1 resonant surfaces in general possess a 
broad spectrum of linearly unstable modes. Moreover, 
the fastest growing mode often has a poloidal mode num- 
ber to ~ O(10). To show this, the q profile shown in 
Fig. ^ is employed, where three q = 1 resonant sur- 
faces are present, at the radii r s i < r S 2 < r s s- By 
evolving the linearized RMHD equations in time, spectra 
of linear growth rates (i.e., dispersion relations) 7n n (wi) 
were obtained as functions of m, as plotted in Fig. for 
Lundquist numbers 5h p = 10 6 , 10 7 and 10 9 . Hereby 
the Prandtl number Pr = SupV has been kept equal to 
unity. Clearly, a variation of S^p (while Pr — 1) retains 
the broadness of the spectrum and 7 max = Max{7r m (TO)} 
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FIG. 5: Evolution of the m = 1 TTM for the q profile given in 
Fig-0 The perturbed kinetic and magnetic energies, denoted 
by _Ek and J5 m are plotted in (a). The growth rates 7k and 
7 m given in (b) are estimated from E k and E m , and 7m ax — 
Max{7n n (m) } is the maximum growth rate in the spectrum in 
Fig- HI for 5h p = 10 6 . The system reaches the fully nonlinear 
saturation in the phase denoted as "annular collapse," i.e., 
t > 200. 

is located at to > 1 in all cases. 

The dependence of the spectrum 7r m (m) on the dis- 
tance between the q — 1 surfaces is most easily inves- 
tigated by considering DTM configurations where two 
q = 1 resonant surfaces are present, located at radii r s i 
and r S 2, a distance D\2 = r S 2 — r s i apart. In Fig. [3] the 
DTM growth rate spectra for profiles with D\2 = 0.05, 
0.1 and 0.2 are shown. While varying D 12 , the local 
magnetic shears at the resonant radii, Si = s(r sl ) and 
s-i = s(r S 2), were not changed. It can be seen that the 
narrower the inter-resonance region becomes, the more 
the poloidal mode number of the fastest growing mode 
shifts to larger values of to, and 7 max increases. On the 
other hand, the growth rate of the to = 1 mode hardly de- 
pends on D\2 and it becomes the fastest growing mode 
for sufficiently large values of D\2. Similar results are 
found for TTMs, whereby it is noted that TTMs tend to 
peak at higher to with higher growth rates than DTMs 
with similar D^. The spectra shown in Figs.|21and[3]con- 
trast with that for the single tearing modes (STM), for 
which the to = 1 mode is dominant and the modes with 
higher to are usually linearly stable. Note that the broad 
spectra of tokamak TTMs that were obtained here are 
similar to those of tearing modes obtained by Dahlburg 
and Karpen for triple current sheets (TCS) in slab 
geometry as a model for adjoining helmet streamers in 
the solar corona. 

It must be emphasized that the dispersion curves in 
Fig. [3 show only the growth rate of the most unstable 
mode for each to. However in general, for a given to 
there are up to three unstable TTM eigenmodes, each 
associated with a resonant surface. To illustrate this, 




FIG. 6: Upper and lower halves of ip (top) and <f> (bottom) 
contour lines in the poloidal cross-section during the annular 
collapse (a) and after the q profile was annularly flattened (b). 

the radial structure of the eigenmodes for m = 1 and 
to = 13 are shown in Fig. which were obtained by 
solving the eigenvalue problem for the linearized Eqs. 
and Here Mi (to) denotes the eigenmode with the 
poloidal mode number to that extends only to the inner- 
most resonant surface r s i. Note that Mi(l) has the same 
mode structure as an STM. Similarly, M 2 (to) and M 3 (to) 
denote the eigenmodes that extend to r S 2 and r S 3 , respec- 
tively. It is also noted that, for to = 13, Mi (not shown) 
is stable, as can be expected from the linear stability of 
STMs with higher to. Similar eigenmode structures are 
found for DTMs, and indeed similar instability charac- 
teristics are also expected for q profiles with more than 
three q = 1 resonant surfaces. 

After perturbing a large number of unstable modes 
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at random poloidal angles, the m = 1 mode evolves as 
shown in Fig. EI The most remarkable feature here is the 
presence of a phase of nonlinearly driven growth. There, 
the energy of the m = 1 perturbed mode grows expo- 
nentially as exp(7drivci)- In this example, 7drivc ~ 0.16, 
i.e., the m = 1 mode grows at a rate which is one 
order of magnitude larger than its linear growth rate 
7iin(^ = 1) = 16 x 1CP 3 . The nonlinear growth rate 
7drivc approximately equals twice the maximum growth 
rate in the spectrum (Fig. EJ), 7 ma x = 0.08, because it 
results from the nonlinear coupling of m and m + 1 mode 
pairs. 

The growth rates shown in Fig. [2] belong exclusively 
to M 3 (m) modes, i.e., the modes extending to the out- 
ermost resonant surface at r = r s3 , since it was assumed 
here that s% < S3 (si = 0.35, s 2 — —0.56, s 3 = 1.20). 
On the other hand, for a q profile with s\ > S3, the 
eigenmode M\(ra = 1) has a higher growth rate than 
Mj,(m = 1). However, also in this case, the highest 
growth rate among all m modes, 7 ma x, is typically several 
times larger than 7ii n (™ = 1) and therefore the nonlin- 
early driven growth of the m = 1 mode still exceeds its 
linear growth. It is concluded that the rapid nonlinear 
growth of TTMs occurs in a wide range of TTM q pro- 
files. Since the sawtooth crash is generally considered 
to be triggered by the onset of an m = 1 mode H, the 
results presented here suggest that the nonlinear growth 
of TTMs (or DTMs for a hollow current profile) is one 
of the possible mechanisms for experimentally observed 
abrupt sawtooth crashes, i.e., the fast trigger. Let us 
note that coupled tearing modes on two nearby resonant 
surfaces can exhibit exponential non-Rutherford growth 
to saturation even for q > 1 0, 0] . 

Finally, the fully nonlinear regime of Fig.[S]is discussed. 
At about t = 200 electromagnetic turbulence starts to de- 
velop in the whole inter-resonance region r s \ < r < r S 3, 
as can be seen in Fig. H3 (a). The fluctuations flatten 
the q profile through simultaneous magnetic reconnec- 
tion in the whole inter-resonance region. This partial 
collapse (i.e., not involving the core) is essentially com- 
pleted around t = 220 ~ 230. A contour plot showing 
the situation after the q profile was annularly flattened is 
given in Fig. 0(b). An important property of an annular 
collapse as the one shown in Fig. [|J] is that no m = 1 
islands needs to form and therefore the displacement of 
the core plasma can be rather small. 

This phenomenon is similar to the experimentally ob- 
served off-axis temperature collapses without precursor 
oscillations. For example, Edwards et al. 0] reported 
JET experiments where rapid sawtooth crashes without 
evident precursor oscillations, but preceded by a partial 
crash in an off-axis region, were observed. 

A possible explanation for such phenomena was pro- 
posed by Buratti et al. in terms of a "purely growing pre- 
cursor" [l^. However, the results in Fig. show that, if 
a q profile with relatively small inter-resonance distances 



D12 and D23 is formed even temporarily, rapidly grow- 
ing nonlinear TTMs can also cause such partial crashes 
without precursors. 

Precursor-free partial collapses during the ramp phase 
of compound sawtooth oscillations were also observed in 
tokamak discharges with hollow q profiles [lH IliL flil \v\ , 
where q = 1 DTMs can be expected. In analogy with the 
TTM case discussed above, our results indicate that non- 
linear growth of q = 1 DTMs can lead to such a partial 
collapse. In cases where the growth rate spectrum peaks 
at lower m (e.g., due to larger Dij or higher Sh p ), we 
have also observed partial collapses (i.e., collapses in the 
inter-resonance region) after which the core displacement 
continues to grow, which is similar to experimental obser- 
vations made in JET, given in Fig. 4 (A)-(C) in Ref. fl2j |. 

In summary, our numerical simulations have shown 
that the simultaneous excitation of unstable 9=1 TTMs 
with high to and their subsequent nonlinear interactions 
lead to a rapid onset of the m = 1 triple tearing mode, 
which qualitatively depicts the fast triggering of sawtooth 
crashes observed in tokamak experiments. Similar behav- 
ior has also been found in simulations of q = 1 DTMs. 
If more than three q = 1 resonant surfaces are formed 
in a tokamak discharge, we also expect similar multiple 
tearing modes that grow rapidly with high poloidal mode 
numbers. We have also presented a scenario, where the 
nonlinear evolution of many unstable TTMs leads to a 
partial collapse of a sawtooth without being preceded by 
an 777 = 1 precursor. Similar phenomena were observed 
during compound sawtooth oscillations in several exper- 
iments |illllElEE3 
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